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In this paper we revisit the Levy-Perdew-Sahni equation. We establish that the relation imphcitly 
contains the conservation of energy density at every point of the system. The separate contributions 
to the total energy density are described in detail, and it is shown that the key difference to standard 
density functional methods is the existence of a general exchange-correlation potential, which does 
not explicitly depend on electron charge. We derive solutions for the hydrogen-like atoms and analyse 
local properties. It is found that these systems are stable due to the existence of a vector potential 
A, related to electron motion, which leads to two general effects: (i) The root of the charge density 
acquires an additional complex phase; and (ii) for single electrons, the vector potential cancels the 
effect of electrostatic repulsions. We determine the density of states of a free electron gas based 
on this model and find that the vectorpotential also accounts for the Pauli exclusion principle. 
Implications of these results for direct methods in density functional theory are discussed. It seems 
that the omission of vector potentials in formulating the kinetic energy density functionals may 
be the main reason that direct methods so far are not generally applicable. Finally, we provide 
an orbital free self-consistent formulation for determining the groundstate charge density in a local 
density approximation. 
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I. INTRODUCTION 

A key innovation in theoretical solid state physics in 
the last fifty years was the reformulation of quantum me- 
chanics in a density formalism, based on the Hohenberg- 
Kohn theorem TJ. Despite initial resistance, in particular 
from quantum chemists, the method has replaced previ- 
ous frameviforks and provides, to date, the most advanced 
theoretical model for the calculation of atoms, solids, and 
liquids. However, its implementation relies on a rather 
cumbersome detour. While the Hohenberg-Kohn theo- 
rem is formulated exclusively in terms of electron densi- 
ties and energy functionals, calculations today are based 
almost exclusively on the specifications given by Kohn 
and Sham one year after the initial theorem was made 
public 0. While this procedure is generally successful, 
and implemented today in numerical methods optimized 
for efhciency (see e.g the ingenious way ionic and electron 
degrees of freedom are treated on much the same footing 
following a method developed by Car and Parinello 3]), 
it is highly inefficient in one crucial conceptual point: If, 
according to the Hohenberg-Kohn theorem, the electron 
density is the only physically relevant variable of the sys- 
tem, then solving the Schrodinger equation, setting up 
the eigenvectors, and computing the density of electrons 
is an operation, which creates a vast amount of redundant 
information. Every information, pertaining to the solu- 
tion of the single-particle Schrodinger equation and the 
summation of single electron charges is discarded at the 
end of every step in the iteration cycle. One could there- 
fore say that more than 90% of the information created 
in today's simulations is actually irrelevant. The ques- 
tion thus arises: Do we have to create this information 
at all, or can we find a more direct way to arrive at the 



groundstate density of electrons without this cumbersome 
detour via the single-particle Schrodinger equation? 

In 1984, Levy, Perdew, and Sahni (LPS) published an 
intriguing relation in which the density of charge is 
described by a second order differential equation without 
any reference to Kohn-Sham orbitals pj. It introduced 
the possibility to formulate the many-electron problem 
with just one quantity, the electron density. In this 
respect it seemed to be more in line with the original 
Hohenberg-Kohn theorem than standard methods in 
density functional theory, based on the Kohn-Sham equa- 
tions [21 • In the last twenty years, developments in this 
field focussed on a search for suitable density functional 
approximations, especially for kinetic energy functionals. 
It is by now an expansive field of research in theoretical 
solid state physics and quantum chemistry. The main re- 
sults of subsequent work on the LPS relation P| are: the 
usage of systematically constructed Harriman orbitals 
in a three-dimensional generalisation the introduction 
and analysis of the Pauli potential 7, 8J; the analysis 
of the uniqueness and asymptotic behaviour of the local 
kinetic energy the properties of the kinetic energy 
density lid : and weighted or averaged density approx- 
imation IT^ Il3j . For a comprehensive review, see 
[T^ . The list is, of course, not complete. In particular 
it does not contain references to the development of the 
Kohn-Sham theory, e.g. the continued improvement of 
exchange-correlation functionals, as this is not the topic 
of the present paper. Here, we want to reexamine the 
original LPS relation and its properties. 

The outline of the paper is as follows: in Section UTI we 
describe the local properties of the LPS equation and find 
that it entails conservation of energy density throughout 
a quantum mechanical system. In Section IIIII we give 
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the results of analytic applications of the LPS relation 
for the hydrogen atom. We find that the effective po- 
tential is the sum of a non-zero Hartree potential and an 
equally non-zero general exchange-correlation potential. 
The effective potential is zero due to cancellation of the 
two contributions. In Section Hvl we analyse the physical 
origin of this cancellation and find that motion of single 
electrons creates a vector potential A, called the inertial 
field due to its relation to electron motion. Due to this 
potential the root of the charge density acquires an addi- 
tional complex phase. The phase-shift is in line with an 
Aharonov-Bohm effect d^. Based on these findings we 
formulate the general problem for an iV-electron system 
in Section El Finally, in Section we discuss the re- 
sults and estimate their importance for the development 
of orbital free density functional methods. 



potentials and energy densities may not be directly acces- 
sible, we derive them in the following from fundamental 
considerations. 

The kinetic energy density t(r) is formulated for inter- 
acting bosonic jjj particles. It can exactly be rewritten 
in a more useful manner, 



t{v) 



<(/7(r),Vp(r),V2p(r)). 

(5) 

The usefulness of this formulation will be seen further 
down. It should be noted that the second term of the 
above expression is the von Weizsacker kinetic energy 
density [l^ ■ For electron charge contained in an infinites- 
imal volume dV around a point r, the contribution to the 
interaction energy between electrons and nuclei E^^ti^) 
will be 



II. LOCAL PROPERTIES 

A closer analysis reveals that the genuine novelty of the 
LPS relation seems to have been disregarded to this date. 
It is the conservation of energy at a local level. This is 
well in advance of the single-particle Schrodinger equa- 
tion or a many-body framework, where energy is con- 
served only globally. 



jp{r)dV, 



(6) 



where M is the number of nuclei and Zi the charge of 
the ith nucleus. The electron-nuclear energy density at 
a point r is then 



£ea;t(r) 



dEext (r) 
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i=l 



\r-R,\ 



Vextir)pir). (7) 



General relations 



We start with the LPS relation. 



rearranging as 

1 V2p(r)i/2 



The effective potential contains the explicit electron- 
electron interaction (Hartree potential, (r)) and a gen- 
eralized exchange-correlation potential (gxc-potential, 
Vgxci^)) which includes the exchange and correlation ef- 
fects, as well as the potential contributions due to Pauli 
p(r)i/2=^p(r)i/2^ (1) exclusion, Veff{r) = VH{r) +Vgxc{'r). The contribution 

to the Hartree energy from the electron distribution be- 
tween two infinitesimal volume dV and dV' around points 
r and r' respectively is 



2 p(r)i/2 
and multiplying by p(r) results in 



+ Vextir) + Veffiv) = ^, 



(2) 



-\p{vf'^V^p{vf'^+v,xt{v)p{v) + v,fiiv)p{v) = pp{v). 

(3) 

We will show that the potential at a point r multiplied by 
the charge density at this point, describes the potential 
energy density. Thus the equation is nothing else but a 
description of energy conservation for every point of the 
system. 



i(r) + eext(r) + £e//(r) = £tot(r)- 



(4) 



Here, we have symbolized the term ppir) by a total en- 
ergy density etot (r) . Each term refers to a corresponding 
energy density. The resulting equation is equivalent to 
the LPS relation, which therefore contains energy conser- 
vation also in its general form. Since the relation between 



d^ Eh {r,r') 



p{r)dVp{r')dV' 



(8) 



Integrating over dV' and rearranging the results leads to 
the Hartree energy density at point r. 

The same behaviour is found for the general exchange- 
correlation energy density, 

Sgxcir) = ^^^^vgxcir)pir), (10) 

although at this point the exact form of Uga;c(r) is not 
known. 

So far we have proved that the left hand sides of Eqs. 
© and Q are the same term by term. In order to prove 
that the total energy density is equal to the chemical po- 
tential times the electron density (etot(r) = pp{^)) we 
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employ the variational principle to minimize the total 
energy. This procedure also allows to determine the rela- 
tionship between the functional derivatives of the energy 
terms, and the corresponding energy densities. Integrat- 
ing Eq. 0) over the whole space the total energy will 
be 



EtoM = T[p\ + E,,t[p] + EH[p\+Eg,,[p] 

M 



Ir-Rd 



d^rvgxcir)pir). 



(11) 



The variational principle, including the condition for the 
total number of electrons N with a Lagrange multiplier 
/i, provides us with the groundstate energy. 



Etot[p]-p[ / d\p{r)-N 



S 

Tp 

This leads to the following result 
SEtot ST SEext , SEh , SEgx^ 



= 0. 
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+ / d^r' ' + Vgxcir) + p{r) 
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It must be noted that the energy term p carries two sep- 
arate meanings: mathematically, it is a Lagrange multi- 
plier; physically, it is also the highest occupied level of 
the groundstate solution of the Kohn-Sham equations as 
well as the negative of the ionization energy in the ex- 
act DFT as was proved earlier by Perdew et al. [l3|- In 
the LPS relation, it describes the eigenvalue of the prob- 
lem, and the change of the total energy, if one electron is 
removed from the system if TV >> 1 y. 

Here, we find the first unconventional implication of 
the LPS equation. Since it is not specified, in the deriva- 
tion of the chemical potential from the LPS relation 0, 
which electron is actually removed from the system, the 
chemical potential must be equal for the removal of any 
single electron. In this case, it is only possible, like in 
standard DFT, to think of a many-electron system as 
composed of a discrete number of electrons with different 
energy levels, if we assume that the system will balance 
the removal of one electron by a change of all electronic 
eigenstates so that the highest electronic eigenstate re- 
mains empty. But from the previous derivation we may 
also conclude that Etot describes the total energy of the 
system. In this case, p must have a double meaning: it 
is, firstly, the energy eigenvalue of the A^-electron system; 
and pN is equal to the total energy. Secondly, it is also 



the chemical potential of the system, or the negative ion- 
ization energy. This result differs from standard DFT, 
where the total energy is a sum containing the discrete 
energy levels of every single eigenstate. 

To understand the difference we revert to the proof by 
LPS that p is the negative ionization energy of a sin- 
gle electron removed from the A^-electron system, if the 
number of electrons is sufSciently large j^]. Suppose now, 
we remove one electron from the TV-electron system. The 
related chemical potential is symbolized hy pN, the neg- 
ative ionization energy of one electron of the A'^-electron 
system. Now, the system contains A'' — 1 electrons. This 
means that the effective potential will be lower than in 
the first case: 



Veff,N-l < Veff,N- 



(14) 



Considering local energy conservation this implies also 
that pn-1 < Mw- 

For calculating ST/Sp we used Eq. @ and the follow- 
ing rule of functional derivatives: 



if T[p] = J d^rt{p{r),Vp{r),V'p{r)), 
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(15) 



Sp dp 9Vp 9V2p' 
Multiplying Eq. by the electron density p(r) gives. 
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(16) 



Comparing this with Eqs. @ , Q , © , Q , © and (jTUl) it 

can be concluded that the energy densities for the differ- 
ent terms are. 
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p(r) = Vgxcir)p{r) 



dVn 



dp 



(19) 
(20) 
0.(21) 



Here, we arrive at an interesting consequence for the 
general exchange-correlation potential. The energy den- 
sity for each separate term is described by the functional 
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derivative times the electron density, even for the most 
compUcated kinetic energy term. If this procedure holds 
also for the general exchange-correlation term, as written 
in Eq. (|21l) , then the general exchange-correlation poten- 
tial does not explicitly depend on the electron density. 
This is in marked contrast to the standard formulations 
in DFT, where exchange-correlations not only depend on 
the density of charge, but are generally parametrized in 
terms of the density (e.g. fl^l. 

The energy conservation on a local level, i.e., at ev- 
ery point r of the system, can thus be directly deduced 
from the LPS relation. Comparing the energy densi- 
ties with the functional derivatives of the correspond- 
ing energy terms has another consequence for the ki- 
netic energy density t(r) and the kinetic energy func- 
tional T[p] = J d^rt{r). Starting from Eq. inte- 
grating over the whole space, performing the functional 
derivative, and multiplying by the electron density leads 
again to the kinetic energy density. From that cycle it 
can be clearly seen that the kinetic energy density can- 
not contain arbitrary terms whose space integral van- 
ishes. This finding is in contrast to previous assumptions 
0,^3) 01- It is justified by the local energy conservation 
of the LPS relation itself. Let us analyse the asymptotic 
behaviour of Eq. (|13|1 as |r| oo. After considerin 



that lim|r|^oo t'e2;t(i") = and lim|r|^oo ^'e//(r) = '4 



we arrive at 
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|r|— >oo op 



lim 



|r|- 



p{r) 



(22) 



which was found by Yang et al. ^ . This means that there 
is no need of approximating the kinetic energy functional, 
it is exact and unique and has the correct asymptotic 
behaviour, although it describes non-interacting bosons. 

The detailed analysis reveals three important proper- 
ties of the LPS equation: 

• The energy density is conserved at a local level; the 
total energy density is therefore constant through- 
out the system. 

• The kinetic energy density is unique, it does not 
contain any arbitrary terms. 

• The general exchange-correlation functional does 
not explicitly depend on the density of charge. 

From a practical point of view we note that self- 
consistency within the LPS framework will be much 
faster to achieve, since the density distribution is much 
more restricted. In particular the requirement that the 
total energy density must be constant at every point of 
the system should make the construction of fast algo- 
rithms quite easy. In this respect it must also be noted 
that solutions of the equation scale with the volume of 
the system, thus the number of atoms: the method is 
therefore a true order-N method. 



III. ATOMIC SYSTEMS 

In the following we analyse the results of the LPS re- 
lation for simple atomic systems containing one electron 
like H, He"*", and Li^+. Previously, we have found that 
the direct consequence of the LPS equation is the energy 
conservation on a local level. 

etot{r) = pp{r) t{r)+Vext{r)p{r)+VH{r)p{r)+Vgxc{r)p{r) 

(23) 

As a next step, we analyse the potentials and energy 
terms which contribute to the total energy. 



Etnf — T + E, 



ext + Eh + Egxc 



(24) 



depending on an assumption about electron density. In 
general, the electron distribution in atoms is spherically 
symmetric, with the nucleus occupying the centre of the 
sphere. Unless otherwise stated we use atomic units 
throughout the rest of the paper. 



A. Potentials and energy terms 

We assume in the following that the electron density 
decays exponentially within the atom p{r) = ae~^^, 
with a,P > 0. This type of electron density is exact 
for hydrogen-like atoms, with the values (hydrogen) of 
a = I/tt,/? = 2. With this ansatz all the relevant vari- 
ables of the atomic system can be explicitly written down. 
The number of electrons is 

d\p{r) = Ana / drr^e''^'' = ^ (25) 
Jo P 

The following quantities are necessary for evaluating the 
kinetic energy density: 



p{r) 

W'p{r) = 



-l3r 



d^p{r) 2 dp{r) 



dr"^ r dr 
Thus the kinetic energy density will be 
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The kinetic energy is 

T = j d^rt(r) = Ana J drr'^er^'' 
a _ iV/32 

The chemical potential is given by 

fi ^ lim — -— = lim — TT ^ TT 

r^oo piV) r^oQ 8 Iv 
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The external potential is Vext{r) = —Z/r, with Z being 
the nuclear charge; the electron- nuclear energy is conse- 
quently 



Eext = J d^rvext{r)p{r) ^ ~4:TraZ J drre~ 
AnaZ (3NZ 



(29) 



The Hartree potential can be obtained using Gauss' the- 
orem via the radial component of the electric field: 
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The Hartree energy is thus 
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The gxc-potential can be obtained from Eq. 1)23(1 : 
t(r) 
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Adding the Hartree potential to the gxc-potential gives 
us all electron-electron related potentials, 



«e//W = VH{r) + Vgxcir) 
We get the gxc-energy as 



2Z- (3 
2r ■ 



(33) 



Egxc = I d^rvgxc{'r)p{r) = Ana / drr'^Vgxc{r)e '''' 
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According to Eq. 124() the total energy will then be 

E.., . ^ _ ^ + A^2^ 
8 2 16 

-f^ (z-^^'-N 
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(34) 



pN. 



(35) 



Comparing the total energy to the kinetic energy in Eq. 
(|?7)l . we find that Etot = T+V = -T. This demonstrates 
the validity of the virial theorem for atomic systems, i.e. 
2T -1-^ = 0. Moreover, it becomes clear that p is an 
average chemical potential, i.e. p — Etot/N, where Etot 
is the sum of the chemical potentials as we build up our 
system from 1 to N electrons: 



N 

Etot = ^ Pi 
i=l 



N 



(36) 



The single electron eigenvalues pi can directly be de- 
termined by DFT. All terms in Eq. (|23|l are now 
parametrized. The atomic systems are described by three 
parameters Z, N, (3, although, it has to be noted that 
the form of the electron density considerably changes for 
> 1 which means the above procedure is only valid for 
atomic systems with one electron. /3 is directly related to 
the total energy (Etot), a fact known from experiments 
and DFT calculations: 



/3 = 



N 



(37) 



Now, we have all the tools to analyse the potentials for 
systems containing one electron {N = 1). 



B. H atom 

The electron distribution of hydrogen can of course 
be determined analytically, using the single particle 
Schrodinger equation. However, to elucidate the phys- 
ical features of the LPS model, it is instructive to choose 
a different route. In case of taking the groundstate (Is^) 
of the hydrogen the three parameters have the following 
values: Z = 1,N = l,Etot = —0.5. The calculated po- 
tentials and densities are shown in Fig. 

The calculated values a = 1/n and /3 = 2 are exactly 
the same as in the analytical treatment. However, we 
gain additional insights into the relation between Hartree 
and exchange-correlation potentials: they are not zero, as 
they would be following the single particle Schrodinger 
equation, but Vgxdr) = —vnir) at every r. This be- 
havior also defines the limit for the effective potential 
'"effir) = Vgxc{r) + vnir) > 0. Physically speaking, it 
indicates that a single electron does not interact with it- 
self, not, as in the conventional picture, because electron- 
electron repulsion (vh) is zero, but because it is cancelled 
by exchange-correlation. For this reason we called Vgxc 
the general exchange- correlation potential, because it is 
inconceivable, within a statistical many-electron treat- 
ment, that it will be non-zero for a single electron. An- 
other important feature of t/ p is that it can be negative. 
However, this feature was already analysed in detail by 
Sim et al. |l3. 
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larger than the repulsive electron-electron interaction. 
The physical cause, invoked here to explain atomic stabil- 
ity, is different from the standard one, where it is thought 
that a decrease of the atomic volume would increase the 
energy uncertainty of the electron and thus increase the 
electron's energy. The standard account is based on the 
Heisenberg uncertainty relations. The reason we arrive at 
a different account within the LPS framework are the lo- 
cal properties. Since within the standard framework only 
global properties are accessible, an argument can only be 
made on the basis of global properties, e.g. volume and 
energy. However, the LPS framework provides direct ac- 
cess to all local properties. Therefore the argument can 
be made on a local basis, via the local interactions of 
electron density and fields of interactions. 

It is tempting to extend the analysis to the general 
case of isolated electrons. In case of = 1 electron if 
Z becomes very small Z ^ 0, the effective potential in 
will remain zero, because 2Z — P = 0. Electron 
charge therefore will become close to homogeneous. But 
also in this case the effect of electron-electron repulsion is 
cancelled because of exchange-correlation. This seems to 
indicate, that regardless of the actual extension of elec- 
tron charge, electrons will not experience any repulsion 
between the charge densities at different regions of the 
electron: electrons are therefore stable entities. 

In the general case > 1 the total energy of an arbi- 
trary system of electrons must be greater than 



FIG. 1; Potentials (top) and energy densities (bottom) within 
the hydrogen atom. Note that Vg^cc is equal, but of opposite 
sign as vh- The effective potential v^ff is therefore zero. 



C. Stability of atoms 



The mathematical condition for Vgxcif) ~ —VH{r) can 
be determined by comparing the terms in Eqs. H30|l and 
(133). It is /? = 2Z. Since, according to Eq. 



NZ' 



(40) 



This value of Etot gives the exact total energy for a sys- 
tem of non- interacting electrons. The total energy in this 
case cannot decrease below Etot because this would cor- 
respond to an attractive interaction between electrons. 

Analysing iJe+ and LP^ ions and the 2s^ excited state 
of the electron in H atom, all systems containing A^ = 1 
electron, Vf,ff = was always found. 



Et. 



Nf3^ 



NZ^ 



Z2 
2 ' 



(38) 



which is always true for A^ = 1. If the total energy Etot 
decreases to E',E' < Etot, then the decay (3, and, conse- 
quently, the effective potential, will be: 



(3 = V^8W> ^-8Et, 



(39) 



Veff{r) = vgxcir) + VH{r) = _ (2Z - /3) < 0. 

The energy Etot is therefore the minimum energy of the 
hydrogen atom; the actual energy is bound from below. 
The atom is thus stabilized by repulsive electron-electron 
interactions, which cannot be more than compensated by 
exchange-correlation interactions. A decrease of the en- 
ergy, and thus an implosion of the atom, is prevented 
by the fact that exchange-correlation would have to be 



D. Summary 

To sum up the results of this analysis of hydrogen-like 
systems: (i) We found in all cases that the Hartree po- 
tential is non-zero. This is in line with the general con- 
cept of DFT, where extended distributions of electron 
charge always have an effect on electrostatic potentials, 
(ii) We also found that the atoms are stable due to the 
peculiar features of the general exchange correlation po- 
tential: it is a negative contribution, equal in magnitude 
to the Hartree potential, which prevents the repulsion 
of electron charge to have any effect. The extent of the 
charge distribution is thus bound from above. But it is 
also bound from below, since the general exchange cor- 
relation potential does not exceed the magnitude of the 
Hartree potential. Hydrogen-like atoms as well as free 
electrons are therefore stable entities. 
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However, at this point the physical nature of this po- 
tential remains mysterious: it does not explicitly depend 
on the charge distribution, as has been generally derived 
from the LPS relation. It can also not be related to ex- 
change correlation in standard DFT, since the systems 
are composed of only one electron. 



IV. THE ORIGIN OF Wg^^ 

At this point we have to look at the physical situa- 
tion from a more fundamental point of view. In general, 
electrons are described by three physical properties: (i) 
Mass, (ii) charge, and (iii) spin. From a fundamental 
point of view, one of these properties must be related to 
the existence and the properties of the general exchange 
correlation potential. It is straightforward to exclude 
charge itself as its origin, since both effects of charge, 
the electron-electron repulsion, and the electron-nuclear 
attraction, are part of the description. Mass, even though 
it fits one of the characteristics of the general exchange 
correlation, i.e. the attractive nature, cannot be responsi- 
ble for it, because the coupling constant for gravitational 
interaction is tens of orders of magnitude smaller than 
the coupling constant for electrostatic interactions. This 
leaves only one viable option: the general exchange cor- 
relation potential must be related to the magnetic prop- 
erties of electrons. 

In general, motion of charge is related to the exis- 
tence of fields. Conversely, fields will affect the motion 
of charge. From this perspective the motion of electron 
charge within a hydrogen atom is likely to create corre- 
sponding fields, which shall have an effect on the energy 
density at a particular point of the system. There is, 
however, one restriction: these fields cannot be equal to 
electromagnetic E and B fields, as this would lead to a 
change of the energy of the atomic system due to radia- 
tion. The field in question can therefore only be a vector 
potential A. The kinetic energy operator used in the 
previous sections does not contain any field components. 
For this reason it most likely gives only a limited account 
of the physical environment. 



A. The modified LPS equation 

The existence of a field A and the ensuing kinetic 
energy operator i (— iV — A)^ will introduce imaginary 
components into the Schrodinger equation. Real func- 
tions (j) ~ p^/^ are thus not sufficient to describe elec- 
trons in this situation. This requires to generalize the 
wavefunction for complex phases and to reinvestigate the 
LPS relation according to this change. The motivation 
for this procedure is to determine the physical origin of 
the general exchange correlation potential. To this end 
we start with a free electron in a three-dimensional po- 
tential well. In this case the electron can be described as 



a plane wave. 



(41) 



with V the volume of the potential well and k the mo- 
mentum of the electron. The above wavefunction results 
in a constant electron density. 



(42) 



From the above it seems to be clear that the LPS equa- 
tion, Eq. ^ can not be complete, since for the free elec- 
tron the square root of the electron density is constant 
which implies that p, has to be zero, which is obviously 
not the case. A similar remedy was su gges ted earlier 
by Wang and Carter, see section VI. 1 of [lj|. Here, we 
extend their results and assume the most general wave- 
function to be complex. The wavefunction contains two 
important features of the electron: amplitude and phase, 
generally 



*(r,t) ^ p(r,t)i/2e''^('"'*). 



(43) 



Writing the Schrodinger equation with this wavefunction 
in the most general way (we omit space and time depen- 
dencies for brevity). 



dip . 1 dp^/^ 
"W^*pi/2 dt 



1 



(44) 



results in two equations for the real and imaginary parts 

dip ivv/' 1 



dt 

1 dp'/^ 
pi/2 dt 



p 

•,1/2 



1/2 



+ l^i^vY + '"ext +VH +'Sie Vg 



VpV2v^ 1 



,1/2 



V -I- 5m V, 



gxc- 



(45) 



Here, the most general Vgxc is a complex potential. In a 
stationary state an eigenvalue equation applies. 



d-^ 



_dip .J_dp^ 
^'dt^^'p^ dt 



= = (46) 



which again results in two equations (real and imaginary 
parts) , 



1 V2pi/2 1 

P = ^1/2 +2^^*^^ + Vext + VH Vgxc, 



P 



1/2 



W (fi + Vg 



(47) 



Eq. H45|) restricts the wavefunction to the simpler form. 



4'(r,t) =p(r) 



1/2 i^p{r) -ifj.t 



ip{r)e' 



(48) 
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The first of Eqs. I|47|l is the modified LPS equation where 
the phase function of the electron explicitly occurs mod- 
ifying the kinetic energy of the electron system (see next 
section). Applying this equation to the free electron gives 
us the correct result for the total energy for TV = 1 
electron) : 



phase 



(49) 



This modification accounts for the well known problem in 
orbital free methods, that the kinetic energy functionals 
of atomic systems, usually described by some modifica- 
tion of the von Weizsacker functional (see above) , or the 
kinetic energy functionals of free electron systems, like in 
many metals, are completely different. Wang and Carter 
proposed a Lindhard function , which connects these 
two extreme cases in a common mathematical model. 
By contrast, we find here that the difference can be ac- 
counted for by the addition of a complex phase to the 
root of the charge density. 

We can sum up the results of this analysis by saying 
that it is essential, for the description of free electrons, 
to take into account both the amplitude and the phase 
information. As a consequence of this, the LPS equation 
has to be modified. The phase information, we found, 
enters as an additional term in the kinetic energy density, 
(V^)V2. 



The physical meaning of general exchange 
correlations 



In previous sections we have shown that it is essen- 
tial, for the stability of atoms, that the general exchange 
correlation is not zero. We claim that this is also the 
case for the free electron, see next section. Furthermore, 
it was shown that a complex phase is essential for ob- 
taining the correct energy eigenvalue. Here, we want to 
determine the origin of this phase and show that the gen- 
eral exchange correlation is due to an inertial field of the 
electron. To this end let us consider the modified mo- 
mentum of a charged particle travelling through a region 
of space with nonzero electromagnetic potentials. It is 
well known that the momentum of such an electron is 



p = -iV - A, 



(50) 



where A is the vectorpotential, while the curl of the vec- 
torpotential is the external magnetic field. 



B = V X A. 



(51) 



The Aharonov-Bohm effect [l^ established the impor- 
tance of electromagnetic potentials in quantum physics. 
For example, in a system containing an external magnetic 
field, e.g. a solenoid, an electron is affected by the vec- 
torpotential, even propagating through a region of space, 
where the external magnetic field is zero. It acquires the 



ipAB = / A(r)dr. 

path 



(52) 



The effect is also observed in experiments |20|, |21[ . Un- 
der the condition that an inertial field A(r), i.e. a field, 
which is not due to external sources, but due to the prop- 
agation of the electron itself, exists, the properties of the 
general exchange-correlation potential can be readily de- 
rived. The Hamiltonian now takes the form, 

H = -A{v)]^ +v,,t{v)+VH{v) (53) 

= -iv2 + lA(r)2 + ^[A(r)V + VA(r)] 

+Wea;t(l") + Vh{y). 

Comparing this to Eq. (^) and using that Vi,ff = vh + 
Vgxc it is straightforward to conclude that the generalized 
exchange-correlation operator is 

vg.civ) = iA(r)2 + ^[A(r)V + VA(r)] 

= iA(r)2 + ^[V.A(r) + 2A(r)V]. (54) 

Using this notation the Hamiltonian can be written as 

H = -iy^ + vext{r) + VH{r) + Vgxc{r). (55) 

The actions of each term of the Hamiltonian on the time- 
space separable wavcfunction, Eq. H48|l are 



Vext{r)'^{r,t) = Uext(r)«'(r,i) 
?)ff(r)^'(r,i) = Wif(r)*(r,t), 
'1 



(56) 



f;g:,c(r)*(r, t) 



A{ry - A(r)Vv3(r) ^{r,t) 



. /V- A(r) Vp(r)i/2 ^ . , , 
-i I + \:,. A(r) *(r,t). 



2 p(r)i/2 

It is important to mention that the kinetic energy is mod- 
ified by a term, [Viy9(r)]^/2 which directly relates to the 
phase of the electron. This differs from the conception 
of e.g. Vignale and Kohn |2^, where the field is related 
to current density. In previous work on orbital free DFT 
this term was never taken into account in the way de- 
scribed here. Adding all terms yields two equations for 
real and imaginary part, respectively, 

A' = ^'ryr + llMr) - V^(r)]^ + vUr) + ^^(r). 



0= (iv 



2 p(r)i/2 
Vp(r) 



1/2 



p(r)i/2 



[A(r)-V(p(r)]. 



(57) 
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The second equation leads to the simple solution 



A(r) = V^(r) 



^(r) - J A(s)ds, (58) 

ro 



where (p{ro) = 0. The phase of the electron is there- 
fore related to a quantity which can be interpreted as 
a vectorpotential. The result seems not surprising, con- 
sidering the Aharonov-Bohm effect 0|. However, an 
objection may be raised that this effect relates external 
A fields to observed phase shifts, and that a phase-shift 
due to propagation of the electron itself will change the 
theoretical predictions for the experiments. The readers 
may convince themselves that this is not the case. If 
two electrons passing on either side of a solenoid possess 
the same inertial field A, then the phase difference be- 
tween any two paths with the same endpoint around the 
solenoid will be given by 

AipAB = j> Aextir)dr + j> A(r)dr = j Be:,t(r)dF -I- 0, 

(59) 

it is therefore determined exclusively by the external 
field. The main difference in our present model is that 
A(r) is the inertial field due to electron propagation, not 
a vectorpotential caused by external fields. The magnetic 
field due to this vectorpotential will vanish, since the curl 
of a gradient is always zero. Substituting Eq. (|58|l into 
the first of Eqs. lf57|) gives 



1 V2p(r) 



1/2 



+ Vext{-r) + Uff(r). 



2 p(r)i/2 

The electron wavefunction takes the form, 



(r,t) = p(r)i/2e'X>(^)%-Mt. 



(60) 



(61) 



Using this solution, we find that the real kinetic energy 
density is not the bosonic kinetic energy density, em- 
ployed in the original LPS relation, but that it possesses 
an additional term, depending on the inertial field A. 

5ietfl(r) = -ip(r)i/2vV(r)^/^ + ip(r)A(r)2,(62) 

^mtniv) = ip(r)VA(r)-f p(r)i/2Vp(r)i/2A(r). 

In this case we obtain a corresponding exchange correla- 
tion energy density described by: 



3fieeg^c(r) = -^p(r)A(r)2, 

3m ega;c(r) = tR{v). 



(63) 



If the inertial field A{y) vanishes, the exchange correla- 
tion potential egxc(r) is zero. In addition, the sum of the 
real kinetic energy density and the exchange correla- 
tion energy density is always equal to the bosonic kinetic 
energy density: 

3?e t{v) + 5Re £g,,(r) - -]^p{vY'^^^ p{vY'^ . (64) 



This entails that the fermi-ionic nature of electrons is 
due to their inertial fields. We shall show further down 
that a gauge transformation, which eliminates the in- 
ertial fields, reverts the problem back to a problem of 
interacting bosons. 



C. Free electrons 

To demonstrate the consequences of the framework de- 
veloped, let us analyse a free electron enclosed in a finite 
volume V . The problem is quite interesting to analyse, 
as DFT does not describe the system correctly. This is 
due to the fact that in DFT both exchange correlation 
potential Vxc and Hartree potential vh vanish, which is at 
odds with the physical situation comprising a finite dis- 
tribution of electron charge. In our description we treat 
the Hartree term in the correct manner. From general 
considerations we have found that Vgxc = ~vh for sys- 
tems of one electron. It should also be noted that v^xt 
is zero within the box and infinity outside. As the phase 
of the wavefunction in Eq. (|41|l is ipir) = kr, the inertial 
field of the free electron is 



A(r) = Vipiy) = k. 



(65) 



This means that the inertial field is just the momentum 
of the free electron system which is uniform within the 
volume. This feature of free electrons was predicted by 
heuristic arguments some years ago 23]. The curl of 
the inertial field is, of course, zero. The magnetic field 
therefore vanishes. Taking this fact and the explicit form 
oivgxc in Eq. I|54|l into account, we can write the modified 
LPS equation as 



-^V2 + 5H(r)+Vc(r) 



ik2+«H(r)-ik^ 



zkr 



V 



(66) 



V 



zkr 



Since Wg// = vu — k^/2 = 0, the chemical potential 
/i — k^/2, and all potentials and also the kinetic energy 
density divided by the density have the same magnitude: 



t 

- =VH 

P 



p. 



(67) 



In addition, they are also equal to the components of the 
total energy of the system, as = 1: 



E 



H 



-Er, 



(68) 



As Eh + Egxc = 0, Etot = T and V = 0. However, 
this is not the general formulation of the problem. In 
the general case, elaborated in detail for a homogeneous 
electron gas in the next section, the field A(r) — Ar(r) 4- 
iA-i (r) will be complex, and we have to account for partial 
waves A"*" and A~ traveling in opposing directions. The 
general treatment will be based on a modification of the 
LPS equation to account for this situation. 
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D. Electrons in hydrogen 



The limit for r — > is consequently: 



Let us determine the inertial field of the electron in a 
hydrogen-like atom. The groundstate density and wave- 
function can be written as 



p{r) 



-2Zr 



TT 



1/2 



(69) 



The system has a total energy of — Acting Vg^c in 
Eq. H54|l on this wavefunction results in two equations for 
the real part and imaginary part, respectively, 



5Re Vn 



\Mr)\ 
V • A(r) 



ZAr{r) = 0. 



(70) 




Here, we assumed that A(r) has only radial components, 
Ar[r). We have already calculated Vgxcir) in Eq. (|22J), 
the parameters in our case are = 1 and (3 — 2Z . Com- 
paring this to the real part of Vg^c, the inertial field can 
be obtained. 



(71) 



The inertial field A in this case is imaginary. However, 
this results is only puzzling, as long as the direction of 
motion remains unconsidered. We have found, for the 
free electron, that A is a vector field, which is parallel 
to the motion vector of the electron. If we assume, in 
case of the hydrogen electron, that the vector of motion 
is radial, then it cannot be unique, because in this case 
the electron distribution cannot be stable. This leads to 
the conclusion, that the standard time-independent solu- 
tion must comprise two separate cases, with the vector of 
motion either parallel, or anti-parallel to the radial vec- 
tor. The field A should therefore be a superposition of 
two partial fields, and A~ . With the ansatz for the 
partial fields. 



yl+(r) = e*x('') A-{r)^-e 



ix{r) 



(72) 



we obtain for the total field Ar{r) 

Ar{r) = A+(r) +A^(r) = 2isinx(r). (73) 

The inertial field in hydrogen is thus described by a radial 
function x(r), which complies with 



/ 1 /2 
x{r) = arcsin I ='=2'\/ 2e^2Zr 



Z 



1 



(74) 



For Z <2 this leads to a well behaved solution, since 
(-2Zr)" 



-2Zr 



l-2Zr + 2Z^r 



2„2 



n! 



(75) 




2e 



-2Zr 



Z 



(76) 



Note that the real components of (r) and (r) pos- 
sess opposite signs, in line with the previous findings for 
free electrons. The real and imaginary components of 
A+(r) are shown in Fig. |21 It is interesting to note 
that A+(r) (i) does not show a singularity like the ex- 
ternal potential, and (ii) does not decay exponentially. 
If we assume that the real part of the inertial field is re- 
lated to electron motion, then the differential of the phase 
V4> = 3?e A+(r) can be used to calculate the phase. Since 
the phase is related to the wavelength of the electron via 
dr/\ — d(i>/2-K, we can also calculate the wavelength of 
the hydrogen electron as a function of radius: 



A(r) 



dr 
d(f){r) 



2tt. 



(77) 



As seen in Fig. |21the wavelength near the hydrogen core 
is about eight or 57r/2 atomic units, it saturates for high 
radii at about 37r/2. From a physical point of view, the 
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FIG. 2: Real (red graph) and imaginary (blue graph) compo- 
nents of the inertial field (r) within hydrogen. The compo- 
nents are equal at the limit of zero radius. The wavelength of 
the hydrogen electron (dashed graph) saturates at large radii. 



real part of A increases as the wavelength decreases. By 
contrast the imaginary part of A decreases with decreas- 
ing density and increasing radius. Tentatively, one could 
thus interpret the real part of A as a physical quantity 
related to motion, and the imaginary part of A as some- 
thing akin to internal friction of the electron. At the 
boundary of the nucleus, or for r — > 0, both components 
are equal, signifying that at this point the energy of mo- 
tion is completely transformed into friction. 
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E. A system of A'^ free electrons: the Pauli principle 

Here, we demonstrate that the inertial field also ac- 
counts for the Pauli exclusion principle in a natural man- 
ner. Separate Pauli potentials, as proposed by March, or 
Levy and Ou-Yang 0, |1| are therefore unnecessary. To 
this end we consider the simplest possible system, the 
system of free electrons. The inertial field of a single 
electron is given by Ai, the corresponding energy of the 
electron is 



1 



1, 



El = -Ai = -k 



(78) 



In a system of two electrons the two inertial fields will 
be superimposed. Assuming, initially, that superimposed 
fields are parallel in space, the combined field of the two 
electrons will be 

A2 = 2Ai E2 = ^2^Aj=4:Ei. (79) 

Since the density of charge will increase by a factor of 
two, the Hartree potential and the general exchange cor- 
relation potential will increase by the same amount: a 
system of two electrons remains thus stable due to the 
superposition of the two inertial fields. A different way 
of expressing this situation would be that for every sin- 
gle electron the effective potential remains zero. For a 
system of N electrons an equivalent description holds, so 
that the total energy of N electrons is given by: 



En = N^Ei. 



(80) 



The result is exactly the same as the one derived for 
a free electron gas in one dimension, which is usually 
accomplished using periodic boundary conditions. Here, 
we have only used the property of any field, i.e. the 
superposition of amplitudes. 

However, in a general three dimensional system the 
direction of fields can vary. In this case the assumption 
of parallel vector potentials is not the state of minimum 
energy. Then the state of minimum energy is described 
by: 



N 



(81) 



Since the individual fields A^ will have equal length 
I Ail = |ki|, the minimum of the total field Ajv will be a 
sphere, equivalent to the Fermi sphere in a three dimen- 
sional free electron model. Note that in all these cases 
the addition of one electron increases the volume of oc- 
cupied states in reciprocal space according to the Pauli 
exclusion principle. The volume of one cell in this sphere 
is determined by the amplitude of the vector potential. 



F. Gauge transformations 

In this section we assume a stationary electron state 
and investigate how it is transformed applying a gauge 



transformation. The effect on the Schrodinger equation is 
also determined. This analysis shall aid us in developing 
methods for an effective solution of the electron problem. 
Here, the aim is to find a suitable gauge for the inertial 
field A. Electromagnetic potentials can be transformed 
with a gauge transformation which results in the same 
electromagnetic fields. 



A'(r) = A(r 

VeLscalir)' = Vel.scalir) + 

E' = E : B' = B. 



vr(r,t), 

dt ' 



(82) 
(83) 



with Vei^scai — Vext+VH the clectromagnetic scalar poten- 
tial and r the gauge function. The wavefunction trans- 
forms like 



^''(r,i) = *(r,t)e-'^(""'*\ 



(84) 



It is easy to show that the Schrodinger equation trans- 
forms according to the following rule. 



(85) 



dt 



dT{v,t) 



dt 



*'(r,t). 



where the initial state ^'(r, t) = ■i/;(r)e~*^* is supposed to 
be stationary. Gauge transformations are potentially use- 
ful in reducing the complexity of the problem and finding 
general solutions. For example, the following transforma- 
tion reverts a general wavefunction back to the root of 
the charge density, 

r(r,i) = f A(s)ds-^t ^ A' = A-Vr = 0, 
^'{r) = pi/2(r), (86) 



which describes a bosonic system as 
1 



Vext +VH - 



(87) 



This is the original LPS equation, see Eq. without 
exchange correlation potentials. It describes interacting 
bosons {veff = vh > 0), and does not comply with the 
Pauli principle. The obvious conclusion from this result 
is that the Pauli principle and the fermionic character of 
electrons is due to their phases. 



V. THE GENERAL PROBLEM 

In previous sections it was established that the gen- 
eral problem of finding the groundstate of an A'^-electron 
system can be described as finding the six components - 
three real, three imaginary - of the complex vector field 
A(r), the distribution of electron density p(r), and the 
chemical potential /z. The complex vector field deter- 
mines the phase of electrons, which is not described by 
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charge density alone. From this perspective the seem- 
ing detour taken in Kohn-Sham DFT is quite under- 
standable. Since charge density alone is insufhcient to 
guarantee phase-continuity throughout a system, it has 
to be modelled by continuous phases of single electron 
states, taking into account electron-electron interactions 
via exchange-correlation functionals. 

Within the present context, this detour is unnecessary, 
if the inertial field can be calculated. That in this case six 
additional values have to be determined for every point 
of the system is not a critical problem. After all, the 
whole computation still scales linearly, even though the 
prefactor might be somewhat larger. In addition, bound- 
ary conditions apply for the inertial field, and the chem- 
ical potential has to be constant throughout the system. 
Given a complex vector field A(r) = Ar{r) + iAi(r), the 
general problem can be formulated in two equations: 



5 [{-v^ 



Ai-Ai + A,V - VA,) 

+Wext +VH - lA P^^"^ = 0, 

A^V - 2A^A0 = 0. 



It is straightforward to separate the different contribu- 
tions and to assign them to potentials within standard 
DFT. The square of the real part of the vector field ac- 
counts for the Pauli principle, as discussed earlier. Thus 
we may define a Pauli potential by: 



vp{v) = \ [Aliv)] 



(89) 



The exchange and correlation potentials are defined in 
standard DFT as the difference of the kinetic energy be- 
tween interacting and non-interacting electrons. They 
correspond in this picture to an operator, vxc{^)i de- 
scribed by: 



standard DFT methods; (iii) a Pauli potential up, which 
depends on the amplitude of the inertial field, in standard 
Kohn-Sham theory it is accounted for by computing the 
solutions for single electron states; and (iv) an exchange- 
correlation potential vxci which depends on the phase 
of the inertial field, described by local densities and their 
gradients in standard methods. Note that within the 
present framework many-body effects are related to the 
inertial field A; in general they are thus mediated by field 
interactions. The effective potential, used in the original 
LPS relation, is thus in fact an operator, given by: 



■i)e//(r) = -Oxc(r) + vp{v) + vh{v). 



A. Homogeneous electron gas 



(93) 



As an example, let us discuss the solution for a homo- 
geneous electron gas. From the auxiliary equation Eq. 
IMl we may infer, as the simplest solution, that the real 
part of the vector potential is zero: 



A,(r) = A+(r)+A,T(r) =0. 



(94) 



From translational symmetry for the external potential 
and the Hartree potential, and from the fact that the 
Hartree potential scales with the density of charge, as 
well as from the translational invariance of the charge 
density p{r) = po, we get: 



Mr) = J d^r'^ 



'3„' 



Po / d^r' 



1 



r — r' 



poa 
(95) 



Thus we arrive at the nonlinear equation for the imagi- 
nary part of the potential, or: 



vxc{r)^l[^A^ (r) + A, (r) V - V A, (r)] . (90) 

The Pauli potential and the exchange-correlation poten- 
tial are not independent from each other, since they have 
to comply with the auxiliary condition: 

[VA,-(r) - A,(r)V - 2A,(r)A,(r)] p'^\r) = 0. (91) 



_[_A2(r)-VA,(r)] -|-apo-Vo = M. (96) 

Choosing again, for translational invariance, a solution 
where A^(r) = Af, we finally arrive at: 



apo - Vq. 



(97) 



Given this condition, the general problem is thus de- 
scribed in a similar fashion to the original LPS equation 
by (the dependency on coordinates is again omitted for 
brevity): 



--V + VXC + VP + Vext + VH - P 



= 0. (92) 



The full potential therefore has four distinct terms: (i) 
an external potential v^xt , which only depends on the dis- 
tribution of positive charge, it is described in every DFT 
framework; (ii) a Hartree potential vh, which depends 
only on the charge distribution, which also is present in 



For a neutral system, where the external potential due to 
the positive background charge and the Hartree potential 
must be equal, or apo — Vq = 0, this amounts to: 



P 



2 ■ 



(98) 



In this case the exchange-correlation potential of the ho- 
mogeneous electron gas is described by —Af/2. From 
a known chemical potential p. the inertial potential can 
therefore be calculated in a straightforward manner. Un- 
der the same conditions the vector field for slowly varying 
densities, which we denote by A/j, will be given by (we 



13 



denote by fj,h the chemical potential of the homogeneous 
electron gas of a specific density p): 



VMh[p(r)] 
|VMh[p(r)]| 



(99) 



Consistent with earlier findings, the exchange-correlation 
potential does not depend explicitly on the density of 
charge. The explicit Pauli potential for the N electrons of 
the homogeneous electron gas can be set to zero, because 
it is implicitly contained in the chemical potential ph- 
Ph[p\- hi turn, can be determined by standard methods 
in DFT. 



B. Local density approximation 

If the vector field is calculated from the homogeneous 
electron gas, then the missing terms in the general prob- 
lem are accounted for. This amounts to a local density 
approximation for the general problem. Assuming a dis- 
tribution of M ions, and a number of N electrons within 
our system, and with A/i(r) from Eq. |^ the general 
problem in the local density approximation is described 
by: 

_ I VV(r) l [Vp(r)]^ 1 
^ ^ 4 p(r) + 8 p(r)2 2 '^^ ^ 

^ [Ah(r)V - VA„(r)] p'^^r) 



M 

E 



z. 



r-R, 



,3 / P(r') 
a r 



r - r' 



(100) 



A;,(r) = ^^2il^[pir)] 
N = I <fr'p{Y'). 



VA*fe[p(i') 
|V/^ft[p(r) 



It is interesting to note that the exchange-correlation po- 
tential contains not only the values of A^(r) and p(r), but 
also their derivatives. This is due to the fact that the in- 
ertial fields are related to the phases of the electrons. 

In this case all terms of the equation depend on the 
local charge and its derivatives. A self-consistency cycle 
then can proceed in the familiar manner by iterating the 
charge density distribution until the chemical potential 
is a minimum. As the only input parameter in the cal- 
culation is the density of charge itself, with the auxiliary 
quantities taken from the results obtained for a homoge- 
neous electron gas, the method is a true order- method. 

However, it is not strictly necessary to remain at this 
level of approximation. Assume, that a groundstate so- 
lution has been found. Then the initial values of p{v), 
Ai(r), and p are known. Subsequently, the solution can 
be refined, by (i) using Eq. |2]to determine the real part 
of the vector potential; (ii) solving the general problem 
with the help of Eq. 1^21 and (iii) iterating the vector 



field and the density distribution until the chemical po- 
tential is a minimum. In this case the solution is the true 
many-body solution of the problem. 



VI. SUMMARY AND DISCUSSION 

Let us briefly summarize the findings in the previous 
sections. It was found that the LPS relation describes lo- 
cal conservation of energy density. It can be generally de- 
rived by variational methods from the total energy func- 
tional, and emerges as the Euler-Lagrange equation of 
the energetic minimum. Its application to hydrogen-like 
atoms showed that these systems are stable, because the 
effective potential, Ug// — + Vg^c, is zero. However, 
the general exchange correlation potential, Wg^c, which is 
equal in magnitude but of opposite sign as the Hartree 
potential, cannot be explained from standard exchange 
and correlation effects for systems containing only one 
electron. An analysis with a more general kinetic op- 
erator including the effect of A fields showed that the 
gxc-potential is due to an inertial field A, related to 
electron motion, which leads to a lowering of total en- 
ergy. The phase of free electron is also due to this in- 

r 

ertial field, Lp{r) = J A{s)ds which has the same origin 

ro 

as the Aharonov-Bohm effect [15j. The inertial field A 
was found to be equal to the momentum k for the free 
electron. A system of N free electrons will be deter- 
mined in reciprocal space by a sphere of finite radius. 
The model in this case includes the Pauli exclusion prin- 
ciple in a natural manner. The general problem can be 
formulated in terms of the local charge density distribu- 
tion and the local complex vector field. From the results 
for a homogeneous electron gas we derived the solution 
for the vector field. This allowed to construct a local 
density approximation for the general problem, which is 
stated exclusively in terms of the density of charge, its 
derivatives, and auxiliary quantities inferred from a ho- 
mogeneous electron gas. 

These results are quite unexpected. Although it is 
probably too early to appreciate the full extent of these 
findings, some general points can be made. 

The existence of inertial fields A, related to electron 
motion, seems to relate to the problem of finding trans- 
ferable kinetic energy functionals in DFT. If, as it seems, 
a crucial part of the physical picture has so far been omit- 
ted, then it is not surprising, that so far every formulation 
has been somewhat less than transferable. It is also im- 
portant to note that the central result of the Hohenberg- 
Kohn theorem, i.e. the proof that the groundstate charge 
density minimizes the total energy, remains untouched by 
our findings. A complex phase added to the root of the 
charge density will not alter the charge density at any 
point. In this case the total energy of a given charge 
distribution is still a minimum, even though additional 
contributions due to the complex phase exist. 

From a practical point of view it could substantially aid 
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in the development of local many-electron frameworks. 
The inertial fields are related to electron motion; as they 
are fields, they should be supcrimposablc and have to 
comply with boundary conditions. Electron charge by it- 
self does not have the same constraints. The well known 
and extensively developed framework of classical electro- 
dynamics may in this case also become useful for the 
study of field properties in quantum mechanical systems. 
The field A is related to magnetic properties of electrons, 
even if, as in case of free electrons or hydrogen-like elec- 
trons, the corresponding magnetic field B is zero. It may 
well be that the ultimate answer to the question, what 
spin actually is, will come from a closer analysis of the 
field under varying external conditions. In addition, the 
finding could lay to rest a century-old problem, that of 
the self-energy of electrons due to electrostatic interac- 
tions. So far, all local electron models could not resolve 
the problem, why electrons would retain their size. Here, 
we found that the effective potential Ve/f = Vgxc+VH = 
for single electrons. This means, that no interactions ex- 
ist between different regions of one extended electron. In 
principle, this allows for any size of electrons. 

From a purely mathematical perspective, and in view 
of direct methods in density functional theory, the proof 
of a successful application to general condensed matter 
problems is still missing. Given that the restrictions this 
relation imposes on the charge density distribution are 
much more stringent than in current models, and that 
any direct method is potentially much more efficient than 
a method based on solutions of the Kohn-Sham equa- 
tions, the results seem promising. Progress, in this re- 
spect, will depend on the development of transferable 
and precise parameterizations of the general exchange- 
correlation potential, or the vector field A. A task, which 



has yet to be accomplished. 

From a physical and more fundamental perspective, 
the emerging picture is quite intriguing. A quan- 
tum mechanical system, which finds its groundstate by 
equilibrating energy density throughout, is in princi- 
ple not very different from a thermodynamic system of 
molecules, which accomplishes exactly the same on the 
level of molecular motion. There is thus no fundamen- 
tal difference, if one goes from the analysis of isolated 
molecules to the analysis of electron charge. Further- 
more, it is quite difficult to explain to a layman, why an 
electron, as a point particle, does not immediately attach 
itself to the nucleus of a hydrogen atom. The present 
framework removes this problem, as well as the seeming 
contradiction, that extended electrons, with their inher- 
ent Coulomb repulsion, are nevertheless stable. 

However, here we should add a note of caution. Even 

though the results obtained so far look very promising, 
the success of the framework will ultimately depend 
on a very pragmatic fact: whether it can be used to 
make accurate predictions of solid state properties, and 
whether these predictions are in line with experimental 
data. 
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